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Abstract. The IceCube Neutrino Observatory has recently placed a new limit on the extremely-high energy 
(EHE) neutrino flux lower than the Greisen-Zatsepin-Kuzmin (GZK) cosmogenic neutrino intensities 
expected in the case of proton-dominated primary composition. This result favors the mixed chemical 
composition models and excludes scenarios with strongly interacting neutrinos proposed so far. Nevertheless 
we argue that it is possible to construct such models consistently with the IceCube data if to extend the 
standard model Lagrangian by adding a leptoquark term. It is notable that the resulting Lagrangian couples 
leptons (not only neutrinos!) to gluons through leptoquark excitations in the subrocess Ig — >■ Iqq so that 
leptons of relevant energies in the Earth as well as in Extensive Air Showers will behave as hadrons. The 
latter may lead, for example, to misinterpretations in identifying the chemical composition of EHE cosmic 
rays. The EHE primaries will seem heavier due to production of secondary strongly interacting leptons in 
agreement with the recent observations. We argue that our model naturally predicts many exotic fenomena 
such as the second knee, centauros, excess of high energy muons (>10 TeV), existence of TeV photons, 
alignment, penetrating cascades. We also argue that a clear manifestation of physics beyond the standard 
model has been observed by the IceCube Neutrino Observatory. 
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1 Introduction 

A really intriguing situation can be observed today in the 
physics of Extremely High Energy Cosmic Rays (EHECR) . 
During many years of observations a great amount of ex- 
perimental information has been collected. A lot of differ- 
ent models have been proposed (for reviews see, for exam- 
ple, [I]). From one hand side the existence of the beautiful 
effect - the Greisen-Zatsepin-Kuzmin (GZK) cutoff [2j[3] 
- is confirmed experimentally, from the other hand there 
are indications that the EHECR chemical composition is 
rather heavier nuclei than protons as recently identified 
by the Pierre Auger Observatory [3]. An opposite picture 
is observed by the HiRes which indicates a proton compo- 
sition |S] . There are also the so-called exotic events whose 
nature is still unknown [51[7]. 

From the existence of the Greisen-Zatsepin-Kuzmin 
(GZK) cutoff one has to unavoidably deduce that there 
should be a bulk of neutrinos (referred to as cosmogenic 
GZK neutrinos) concentrated around the threshold of the 
following reaction 



P7 -> A'^ 



(1) 



In the laboratory frame the threshold is equal about 
5x1019 eV. 

It was also realized that cosmic rays beyond the GZK 
cutoff may also exist if they are initiated, for example, by 
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neutrinos (lOj . Models based on the idea (primary neu- 
trinos) have been proposed. For example the so-called 
strongly interacting neutrino models describe the experi- 
mental results very convincingly and seem to be rational 
provided the cross section of neutrino interactions with 
matter grows up to values of the order of 1 mb at such ex- 
tremely high energies ( 10"'^* eV and above). This is much 
higher than the standard model predictions. 

The IceCube Neutrino Observatory has recently placed 
a new limit on the extremely-high energy (EHE) neu- 
trino flux lower than the GZK cosmogenic neutrino in- 
tensities expected in the case of proton-dominated EHE 
cosmic rays [TTlll2j . This result excludes models beyound 
the standard model proposed so far. 



2 Involving leptons in strong interactions 

We stress that the following two states are taken as prin- 
ciples: 

(i) The GZK cutoff exists for protons. 

(ii) The primary EHECR arc protons. 

Models with strongly interacting neutrinos provide a 
rational explanation of the shape of the cosmic ray spec- 
trum above 10^^ eV. They suppose that the GZK cos- 
mogenic neutrinos initiate extensive air showers, just as 
hadrons and the observed spectrum is nothing but a su- 
perposition of two components initiated by the primary 
GZK neutrinos and protons respectively (see, for example. 
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Fig. 1. Lepton-gluon interaction through leptoquark excita- 
tion (S - skalar and V - vector leptoquarks) 



[13] )• AH such models proposed so far have been convic- 
ingly excluded by the IceCube Neutrino Observatory. 

How to consistently involve leptons in strong interac- 
tions? We propose here a method using a leptoquark La- 
grangian. The existence of leptoquarks was predicted by 
W. Buchmuller, R. Riickl and D. Wyler in 1986 nHI. 

Let us extend the standard model Lagrangian [Lsm] 
by adding a leptoquark term (Clq) 



^ — ^SM + ^LQ ■ 



(2) 



The general SU{3) x SU{2) xU{l) invariant Lagrangian 
that couples neutrinos to leptoquarks has the form (14) : 



^LQ = gi qliT2hSl + gz q'iiT2TlL^3 

+92 d''Rr^LV2^, + 52 u^7^^lV^2m + h2 rfi?/Li?2 

+ h2 URiLR2 + hi QLj^iLUi^ + h3 qLT-f^ (.L^^t, 

+ h.c. . (3) 

The Lagrangian ^ contains eight scalar and vector 
leptoquarks per family grouped into weak isospin doublets 
or triplets (indicated by indexes 1 - singlet, 2 - doublet, 3 - 
triplet) and couples leptons to gluons through leptoquark 
excitations as shown in Fig. 1. 

Therefore one can calculate the corresponding cross 
sections of the following reactions: 



Ig qL, 



(4) 



where I are leptons {vi, l^), g - gluons, L - leptoquarks. 
The leptoquarks subsequently decay into a lepton-quark 
pairs so we have the following subprocesses: 



kg qL ^ qqlj, 



(5) 



The advantage of this model in comparison with other 
ones proposed so far is the fact that not only neutrinos 
but also charged leptons interact strongly at the appropri- 
ate energies, probably near the leptoquark pole (resonance 
production) . 

The final step is to convolute the cross sections with 
gluon densities of the nucleon. 

It is enough to show that within this model one can 
obtaian lepton-nucleon cross sections of the order of 10 
mb to be in agreement with observations of EHECR. 

To demonstrate how this model works, let us consider 
the excitation of the vector leptoquark from the singlet, 
i.e., consider the effective Lagrangian: 



hiqLinLU, 



(6) 



In the leading order we have obtained the following 
cross section: 



a(s) 



fiiy) + .f2{y)\n 



(7) 



where 



/i(y) = 2(l-i)(l + A-llny), 

y y y 



f2{y) 



1 



Ml and are masses of the Icptoquak and the quark, 
as ~ the strong coupling constant, y = s/Mf^, s - the 
center of mass energy, A - a constant corresponding to / — 
L — q vertex. Our model contains two parameters beyond 
the standard model: A and M^. 

Now we convolute the cross section ([7]) with the gluon 
distribution of the nucleon: 



a{E) 



a{xs)g{x, s)dx, 



(8) 



where E is the incident lepton energy in the laboratory 
frame. We choose a usual parameterization of the gluon 
distribution: 



(1 



(9) 



We have taken m=5, n=1.5, ag = 0.1, A=l, M=500 
GeV and performed calculations for two values of m^: 
uiq = 0.2 GeV and = 1 GeV, xq = M'^/As, A - is 
the mass number of the air/earth nucleus which contains 
the target nucleon (see Fig. 2). 

It should be noted that the gluon density at such ex- 
tremely high energy transitions is unknown. One can sup- 
pose, for example, that the leptoquark production in lepton- 
nucleus interactions (in the air or in the Earth) becomes 
possible due to the gluon density is very high at such a 
regime. It is also very important that the incident particle 
is just a lepton for reaching the region of ultra low values 
of the Bjorken variable. In pp collisions it would be more 
difficult to probe the same values of x since the condition 
AP « X1X2S takes place. 

We emphasize that our numerical calculations must be 
regarded only as a qualitative demonstration of the model. 
There are a number of unknown parameters, apart from 
the leptoquark masses and the constants A. 

Nevertheless the model proposed in this paper has an 
impressive predictive power for astrophysics. We demon- 
strate this fact in the next section. 



3 Predictions 

It is notable that one can make a number of significant 
predictions without numerical calculations. 
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Fig. 2. The leading order cross section of a lepton-nucleon in- 
teraction through vector leptoquark excitation v.s. lepton en- 
ergy in the laboratory frame. Two values of quark masses are 
considered: mq=0.2 GeV (blue curve) and 1 GeV (red curve), 
M is taken to be 500 GeV 

1) The recent results of the IceCube Neutrino Observa- 
tory favor the idea that the EHE leptons interact strongly 
and therefore may be absorbed in the Earth. 

2) Within the proposed model, it is possible to explain 
the energy spectrum of EHECR as a superposition of two 
components initiated by the primary GZK neutrinos and 
protons respectively. 

3) The EHE primaries will seem heavier due to pro- 
duction of secondary strongly interacting leptons. For ex- 
ample some muons may be misidentified as hadrons. 

4) The TeV photons may be explained as bremsstrahlung 
of EHE charged leptons (say electrons) produced in the 
first interactions high in the atmosphere. 

5) Centauros, excess of high energy muons (>10 TeV), 
penetrating cascades seem not to be so exotic in the model 
of strongly interacting leptons and can be easily explained. 
For example, the probality that the neutrinos with ener- 
gies out of the leptoquark pole can penetrate deeply in the 
atmosphere and then interact strongly is not equal to zero. 
Centauros seem to be nothing but leptons misidentified as 
hadrons. 

6) The alignment may exist due to leptoquark decays 
L ^ ql (only one hadronizing quark in the final state) . 

7) The events observed above the GZK cutoff and de- 
bated so far may be interprited as induced by neutrino 
primaries. 

8) The IceCube data favor the assumption that all the 
three families of leptoquarks exist since no EHE tau lepton 
has been observed. 

9) The masses of leptoquarks are of the order of 1 TeV. 

In conclusion we would like to note that all the ar- 
guments given in this section suggest that we deal with 
physics beyond the standard model and its clear manifes- 
tation has been observed by the IceCube Neutrino Obser- 
vatory. 
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